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ABSTRACT 

We report multi-epoch GMRT HI observations of the z = 0.3127 damped absorber 
towards the quasar PKS 1127-145 , which reveal variability in both the absorption 
profile and the flux of the background source, over a time-scale of a few days. 

The observed variations cannot be explained by simple inter-stellar scintillation 
(ISS) models where there are only one or two scintillating components and all of 
the ISS occurs in the Galaxy. More complicated models where there are either more 
scintillating components or some of the ISS occurs in the ISM of the z = 0.3127 
absorber may be acceptable. However, the variability can probably best be explained in 
models incorporating motion (on sub-VLBI scales) of a component of the background 
continuum source, with or without some ISS. 

All models for producing the variable 21cm absorption profile require small scale 
variations in the 21cm optical depth of the absorber. The length scale for the opacity 
variations is ^ 0.1 pc in pure super-luminal motion models, and ^ 10 pc in pure 
ISS models. Models involving sub-luminal motion, combined with scintillation of the 
moving component, require opacity variations on far smaller scales, ~ 10 - 100 AU. 
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1 INTRODUCTION 

The flat spectrum radio- loud quasar PKS 1127-145 
{z ~ 1.187) is known to vary at many wave bands and 
over widely differing time scales. Long term flux monitoring 
(e.g. Wehrle et al. 1992; Bondi et al. 1994) has shown that its 
flux at several radio frequencies varies on time scales rang- 
ing from days to months. The source is an intra-day optical 
variable ( Romero, Cellone fc Combi 1999 ). It has also been 
detected in the X-ray (^ic j3ert_et^ al. 199S ) and the gamma- 
ray bands (Thompson et al. 1995). At energies > lOOMeV, 
EGRET observations found flux variations between the two 



of the core at the third epoch, suggestive of the emergence 
of a new component. 

Damped Lyman-a absorption was detected at z ~ 
0.3127 towards PKS 1127-145 by Rao & Turnshek (1999), 
using the Hubble Space Telescope. The column density of 
the absorbing system was found to be quite high, Nhi = 
5.1 X 10 per cm . Optical observations have so far been 
unable to provide an unambiguous identification of the ab- 
sorber, as there are three galaxies close to the line-of-sight, 
at redshifts very similar to that of the d amped absorption 
(Lane et al. 1998; Bergeron & Boisse 1991). 21cm absorption 



epochs of observation (McLaughlin et al. 1996). 

Sub-milliarcsecond im aging of PKS 1127-145 at 15 GHz 
( Kellermann et al. 199^ ) has shown that it has a simple 



core-jet structure, with a prominent knot at the location of 
a sharp bend in the jet. This core-jet stru cture is also seen 
in VLBI observations at lower frequ encies (Tingay, Murphy 
fc Edwards 1994 iBondi et al. 1996|). Multi-epoch VLBI ob- 



was also detected from this z = 0.3127 system by Lane et 
al. (1998), using the Westerbork Synthesis Radio Telescope 
(WSRT). Higher sensitivity Giant Metrewave Radio Tele - 
scope (GMRT) observations (Chengalur fc Kanekar 200C) 



servations at 1.67 GHz (pondi et al. 1996) showed a strong 



increase in both the flux density as well as the compactness 
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showed that the absorption profile has a number of nar- 
row components, spread out over a velocity range of ^ 100 
km s~^. 

We report here on multi-epoch GMRT observations of 
this 21cm profile. As described in section ^ below, the 21cm 
absorption profile is observed to vary on a timescale of a 
few days. This is the second extra-galactic 21cm absorber 
to show variability, the other system being the z ~ 0.524 
damped absorber towards the BL Lac object AO 02354-164 
(Wolfe, Briggs fc Jauncey 1982). We discuss, in section ^ 
models which might produce the observed variations. 
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2 OBSERVATIONS AND RESULTS 

PKS 1127-145 was observed with the GMRT at 8 epochs 
during April - May and September - October 1999. The 
backend used was the 30 station FX correlator, configured 
to give 128 channels over a total bandwidth of 1 MHz, i.e. a 
channel spacing of ~ 7.8 kHz (or - 2 km s"^). The number 
of antennas used for the observations varied between 7 and 
12. The varying baseline coverage is, however, of no conse- 
quence since PKS 1127-145 is an extremely compact source 
(angular size Si 15 milliarcseconds) and is not resolved by 
even the longest baselines of the GMRT. The observing de- 
tails are given in Table ^. 

Flux and bandpass calibration were carried out using 
3C286, which was observed at least every 40 minutes dur- 
ing each observing session. The data was converted from 
the telescope format to FITS and analyzed in AIPS using 
standard procedures. The task UVLIN was used to subtract 
the continuum emission of the background quasar; maps 
were then produced and spectra extracted from the three- 
dimensional data cube. Spectra from different days were cor- 
rected to the heliocentric frame, scaled to the same flux level, 
and then combined to produce the final averaged spectrum. 

The lower panel of Figure |l| shows the final spectrum, 
averaged over the 8 epochs of observation; here, flux is plot- 
ted against heliocentric velocity, centered at z = 0.3127. All 
spectra were scaled to a common continuum flux of 6.3 Jy, 
before they were combined together, i.e. the y-axis is effec- 
tively the optical depth. The RMS noise level on the spec- 
trum (note that this is computed from absorption-free re- 
gions in the final averaged spectrum) is ~ 3.5 mjy. The 
upper panel of the figure shows the RMS across the eight 
epochs of observation, plotted as a function of heliocentric 
velocity. It can be seen that the noise levels are far higher 
at line locations than in absorption-free regions, as would 
be expected for a time- variable absorption profile. Figure ^ 
shows a plot of difference spectra, obtained by subtracting 
the final averaged spectrum from the spectra of individual 
epochs (each scaled to a common continuum flux of 6.3 Jy; 
the y-axes are thus effectively the differential optical depth). 
The average spectrum is also plotted on the two lowest pan- 
els, for comparison. Note that the scale on the difference 
spectra has been expanded by a factor of 10 as compared to 
that on the average spectrum. Finally, the 3a noise on each 
difference spectrum (computed in absorption-free regions) is 
indicated on the left, by the vertical barred lines. 

An examination of figure |^ shows that the variation of 
different components is not synchronized, i.e. the observed 
variability is not due to some simple problem in the scaling 
of the continuum flux. There is no evidence for interference 
corrupting the spectra. Since the spectra show a measurable 
Doppler shift over the period of observation, an interference 
signal, which would stay at a flxed frequency, would be easy 
to pick out. At all our observing epochs, we have observed 
both polarizations; the spectra from the two polarizations 
(which are processed through largely independent receiver 
chains) are identical to within the noise. On one epoch (10th 
May 1999), we have two contiguous observations with the 
line profile shifted by a non-integral number of channels; 
the profiles obtained from these two observations are again 
identical within the noise. Finally, we have verified through 



simulations that the observed variations are not the result 
of finite channel resolution. 

Fitting of multiple Gaussians to the average spectrum 
and to the spectra of individual epochs shows that the 
changes in the profile are consistent with variations in only 
the depths of the various components, with no changes in 
their positions or widths. A total of nine Gaussians were fit 
to the final spectrum. Although the changes in the depths 
of some of the components were found to be correlated (no- 
tably the three deepest components in Figure 0), not all 
the components vary in a correlated manner. The observed 
variations in the 21cm profile of the z — 0.524 damped ab- 
sorber towards AO 0235-1-164 ( |Wolfe et al. 1982| ) can also be 
similarly explained by changes solely in the amplitudes of in- 
dividual components of the profile (i.e. with no changes in 
their centroids or widths). Note that although those authors 
quote a characteristic time scale of six months for changes 
in the z — 0.524 absorber, examination of Figure (5) in 
Wolfe et al. (1982) shows statistically significant differences 
between spectra taken as little as 2 days apart. 

Table ^ also lists the measured flux values of PKS 1127- 
145 for each observing session, along with error bars. The 
fluxes are found to vary widely, with values ranging from 
~ 4.0 Jy to ~ 7.9 Jy. The most dramatic change is between 
the 8th and the 10th of October, a flux increase of 1.8 Jy over 
two days. Our experience with the GMRT indicates that the 
flux calibration is reliable to ~ 15%, in this observing mode. 



3 DISCUSSION 

As described above, the observed 21cm optical depth (r2i) 
of the z = 0.3127 absorber varies over timescales of a few 
days. A change in the observed optical depth could arise 
from a number of reasons : (1) changes in the physical con- 
ditions (i.e. spin temperature, column density) in the ab- 
sorber, (2) inter-stellar scintillation of one or more compo- 
nents of the background source and (3) motion of a source 
component across the sky so that the line-of-sight to this 
component traces regions of different opacity. We will, in 
this section, discuss each of the above possibilities in greater 
detail. Models for producing the variable 21cm absorption 
towards AO 0235-1-164 are discussed by Briggs (1983). 

The short timescales involved make it very unlikely that 
the observed variability originates in changes in the column 
density along the line-of-sight, arising from motion of the ab- 
sorbing clouds. Even if one assumes that the characteristic 
scale on which the opacity varies is 20 AU (see below) and 
that the absorption occurs against a highly compact back- 
ground source, variations over a timescale of 2 days would 
require the absorbing material to have transverse velocities 
of ^ 0.1c. This is quite unrealistic for cold gas in a galaxy or 
for tidal debris. Similarly, there is no plausible mechanism 
for recurrent changes of the spin temperature of a cloud of 
this size over a timescale of days. We hence do not consider 
this model further. 

Changes in the line profile could be explained if some 
part of the background continuum flux arises from a region 
of such small angular size that it gives rise to inter-stellar 
3cintUlation (ISS) ([Pennett-Thorpe fc dc Bruyn 2000|; [Rick' 
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Table 1. Observing details. 



Observation epochs 


Centre frequency 


Flux density" 


(1999) 


(MHz) 


(Jy) 


27 Apr 


1082.000 


5.2 ±0.8 


6 May 


1082.000 


7.9 ±1.2 


10 May 


1082.049 


6.1 ±0.9 


10 May 


1081.949 


6.0 ±0.9 


17 May 


1081.949 


7.0 ±0.4 


22 Sep 


1082.050 


4.0 ±0.6 


8 Oct 


1082.050 


4.3 ±0.6 


10 Oct 


1082.050 


6.1 ±0.9 


15 Oct 


1082.050 


7.4 ±1.1 



"The errors quoted include systematic errors, based on earlier GMRT observations of standard calibrators in this observing mode. 
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Figure 1. The lower panel shows the final averaged spectrum to- 
wards PKS 1127-145 plotted versus heliocentric velocity, centred 
a,t z = 0.3127. The RMS noise level is plotted in the upper panel. 



ett 1990; Kedziora-Chudczer et al. 1997 ) on passing through 
the interstellar medium of the Milky Way. The spectrum as- 
sociated with this region will vary, while that associated with 
the non-scintillating flux will not; this will result in changes 
in the net optical depth. ISS will also produce variations 
in the continuum flux; further, these variations would be 
related to the changes in the optical depth. Since we mea- 
sure bo th the background flux and T21 , one can test whether 



Figure 2. Difference spectra on the 8 epochs of observation, plot- 
ted as a function of heliocentric velocity, centred at 2 = 0.3127. 
The panels are labelled by the dates of observation. The lowest 
two panels contain the final averaged spectrum, for comparison. 
Note that the scale on the y-axis is expanded for the difference 
spectra. The 3cr noise on each difference spectrum (computed in 
absorption-free regions) is indicated on the left, by the vertical 
barred lines. 



cussed in more detail below, the lack of correlation between 
the variations in flux and opacity rule out the simplest such 
models. 



15 GHz VLB I observations of PKS 1127-145 (Keller- 



the obssrvations are consistent with ISS models. As is dis- 



mann et al. 1998) have shown that it has a simple core-jet 
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structure. The core is unresolved at sub-milliarcsecond res- 
olution. The simplest scintillation model is that of a single 
compact scintillating component in the core, with the rest 
of the flux not giving rise to scintillation. In this case, how- 
ever, two epochs with the same observed continuum flux 
should have the same observed spectrum; further, the line 
depths should be deeper at epochs of higher flux and vice 
versa. Table (1) shows that the flux values agree within the 
error bars on the 22nd of September and the 8th of Octo- 
ber; however, the spectra at these epochs are considerably 
different. Similarly, the background flux is the same on the 
10th of May and the 10th of October while the spectra are 
again quite different. Given our error bars, the flux change 
required to produce the observed differences in the 21cm 
spectra would have been easily detectable, unless the line- 
of-sight to the scintillating component also coincidentally 
happened to have an optical depth of order unity. The lat- 
ter situation is also, however, not viable since the epoch with 
the highest background flux (5* = 7.9 ± 1.2 Jy, on the 6th 
of May) has the weakest absorption. Thus, the simplest ISS 
model appears to be ruled out. 

The next simplest ISS model is one in which there are 
two scintillating components and a steady component. The 
second scintillating component could, for example, arise in 
the bright unresolved knot in the jet of PKS 1127-145 , 
which is separated from the core by ~ 2.5 milliarcseconds. 
If the optical depth towards each of the scintillating com- 
ponents is comparable to the average optical depth over the 
entire source (i.e. T2\ ~ 0.1) then, as discussed below, this 
model is also ruled out by the data of 22nd September and 
8th October. The difference spectrum on 8th October con- 
sists of two components which are ~ 50 mjy deeper than the 
average spectrum and three which are ~ 30 mJy deeper than 
the average. On the 22nd of September the first two compo- 
nents are weaker (by ~ 30 — 40 mJy) than the average, while 
the remaining three agree (within the noise) with the aver- 
age spectrum. If one assumes that the absorbing clouds are 
distributed such that the changes in the first two line compo- 
nents arise from one of the scintillating source components 
(total change in line depth, AS' ~ 90 — 100 mJy, per compo- 
nent) and the changes in the remaining three stem from the 
other (total change in line depth, AS ~ 30 — 40 mJy, per 
component), the changes in line depth would correspond to 
a flux change of at least 1.3 Jy (for T2\ ^ 0.1), which is not 
seen (other distributions of the clouds would exacerbate the 
problem). Since these two compact components contain far 
more than 10% of the total flux, it is rather unlikely that 
the optical depth in front of the scintillating components is 
~ 1, while the average optical depth is ~ 0.1. 

One could construct more complicated models, wherein 
each of the compact VLBI components contain more than 
one scintillating component. The relation between the ob- 
served spectrum and the measured flux becomes more and 
more complex with an increasing number of scintillating 
components and it would hence be possible to contrive situ- 
ations where both the variations might be explained purely 
as a result of ISS in the Galaxy. However, from the Taylor- 
Cordes (1993) model of the distribution of ionized gas in 
the Galaxy, the scatter broadening expected towards PKS 
1127-145 is 0.4 milliarcseconds. It thus seems unlikely that 



uncorrelated ISS would result, even if the core or knot has 
multiple compact components. We note, however, that ISS 
might also occur in the galaxy (ies) giving rise to the absorp- 
tion at z — 0.3127, while all the above models only consider 
the case of it occurring in the Galaxy. Scintillation could 
prove to be a viable model if significant ISS does occur in 
the interstellar medium of the absorber. 

Motion of a source component across the sky might 
explain the changes in absorption profile, if the absorbing 
system has opacity variations transverse to the line-of-sight. 
The line-of-sight to the moving component samples different 
paths through the absorber at different epochs, resulting in 
changes in the profile. Super-luminal motion on sub- VLBI 
scales has oft en been invoked to explai n rapid variability of 
QSO fiuxes (Rees 1966; Woltjer 1966). If these variations 



are intrinsic, they are difficult to explain in the standard 
context of incoherent synchrotron emission as causality ar- 
guments would imply brightness temperatures far in excess 
of the limit of 10^^ K set by the inverse Compton catastro- 
phe. Super-luminal motion of a source component produces 
a Doppler boosting of the observed brightness temperature 
by a factor over the true source brightness temperature, 
where S is the Doppler factor. Thus, if S is sufficiently large, 
the true brightness temperature would be lower than the 
inverse Compton limit. In these models, the timescale of 
variability can be used to place a lower limit on 5. Doppler 
factors of order 10 are sufficient to account for the vari- 
ability seen in the majority of compact sources; similar val- 
ues have also been measured from VLBI observations of 
super-luminal motion. However, observed variations in the 
1420 MH z flux of AO 0235+ 164 on timescales of days yield 
S ^ 100 (Krauss et al. 1999). Similarly, long term monitor- 
ing of the flux of PKS 1127-145 at 4.8, 8.0 and 14.5 GHz 
(under the Michigan Monitoring program) have found short 
timescale variability which, if assumed to be intrinsic, would 
require 5 ~ 100. Our own continuum measurements (albeit 
with larger error bars) also require S values of the same or- 
der. 

As mentioned earlier, the change in absorption profile 
in this model also requires the absorbing system to have 
opacity variations transverse to the line-of-sight. The extent 
of inhomogeneity in the absorber can be estimated from the 
fact that the physical distance I travelled by the moving 
component between two observations separated by an inter- 
val Atobs (in the observer's frame) is ~ ScAtobs, where we 
have assumed 5 « 7, the Lorentz factor. The correspond- 
ing projected distance in the absorber is ScAtotsiDc/ Ds), 
where Ds and Dc are the angular diameter distances of 
the source and the absorbing clouds respectively. Variations 
in the line profile over Atobs ~ 2 days would require the 
absorbers to have opacity variations on similar scales, i.e. 
~ 1.6 X 10~^S{Dc/Ds) parsecs. Using the source and ab- 
sorber redshifts in a flat, matter-dominated FRW cosmol- 
ogy (with Ho = 75 km s~^ Mpc"^) yields Dc ~ 775 Mpc 
and Ds '-^ 1184.5 Mpc. The absorbers must thus be inho- 
mogenous on length scales of ~ 0.1 parsecs to explain the 
changes in the line profile in the above model. (Note, how- 
ever, that the coherence scale of velocities in the absorber 
must be much larger, to account for the unchanged posi- 
tions of the velocity centroids of individual line components.) 
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The Doppler factors required in this scenario are, of course 
rath er large but might be explicable in shock-in-jet mod- 



els ( Marscher fc Gear 1985 ; Qian et al. 1991 ; Blandford & 
Konigl 1979p , which are variants of the standard relativistic 
beaming paradigm. 

Sub-luminal motion of a source component could also 
explain the variations in both the flux and line proflle, if 
the moving component was sufBciently small so as to give 
rise to interstellar scintillation. Of course, there would be no 
correlation between the measured flux and the line profile 
in this situation, since the line-of-sight to the scintillating 
component would itself change. This scenario would, how- 
ever, require the absorbing clouds to have structure on very 
small scales. For example, speeds of the order of 0.1 c (at 
z = 1.187) would imply changes in opacity on scales of ~ 
20 AU, at the redshift of the absorber, to account for the 
changes in line profile over two days. Such fine structure in 
HI opacity, on scales of 10 — 100 AU, has, in fact, been seen 



earlier in our Galaxy, in VLBI 2 1cm absorption studies (Fai- 
son et cl. 1995; Frail et al. 1994). In fact, these observations 
have found the optical depth in the Milky Way to vary by 
a factor of two, over the above length scales; this is more 
than sufficient to explain the changes seen in the absorption 
profile of the z — 0.3127 absorber. We note in passing that 
several spin temperature measurements, in particular those 
of damped Lyman-a systems would become suspect, should 
10 - 100 AU be the typical size scale of HI "clouds" . As em- 
phasized by Deshpande (2000), however, opacity variations 
on scales of 10 ~ 100 AU do not neccassrily imply the ex- 
istence of physically distinct extremely dense clouds of this 
size. If these variations are only random opacity fiuctuations 
on small scales, with a much weaker systematic gradient in 
opacity, the spin temperature measurements are probably 
reliable. 

The final model we mention is one where gravitational 
microlensing of a moving background component by stars or 
stellar clusters in the z = 0.3127 absorber gives rise to varia- 
tions i n both the fi ux and the line proflle ( Gopal-Krishna & 
Subran anian 1991). Microlensing causes the apparent trans- 
verse velocity of a moving component to increase; this would, 
in turn, increase the length scale on which opacity variations 
are required to exist in the absorbing source. Microlensing 
scenarios are strongly dependent on the geometry and struc- 
ture of both the source and the lens. Further, a source which 
is compact enough for micro-lensing to be important would 
probably also give rise to ISS; we will hence not explore any 
detailed microlensing models here. We note, however, that 
microlensing is a likely scenario, given the fact that an ab- 
sorbing galaxy lies along the line-of-sight to the background 
quasar; simultaneous multi- frequency observations could be 
used to test the presence of microlensing, due to the achro- 
matic nature of this phenomenon. 

In summary, while there are a number of models which 
could produce the observed variability in both the flux and 
the 21cm profile, all models require the optical depth of the 
absorber to have small scale structure. The length scale for 
the opacity variations is ~ 0.1 pc in pure super-luminal mo- 
tion models, and 10 pc in pure ISS models. Models involv- 
ing sub-luminal motion, combined with scintillation of the 



moving component, require opacity variations on far smaller 
scales, ~ 10 - 100 AU. 
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